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ABSTRACT 

Hot subdwarf B stars (sdBs) are evolved core helium-burning stars with very 
thin hydrogen envelopes. In order to form an sdB, the progenitor has to lose 
almost all of its hydrogen envelope right at the tip of the red giant branch. In 
binary systems, mass transfer to the companion provides the extraordinary mass 
loss required for their formation. However, apparently single sdBs exist as well 
and their formation is unclear since decades. The merger of helium white dwarfs 
leading to an ignition of core helium-burning or the merger of a helium core and 
a low mass star during the common envelope phase have been proposed. Here 
we report the discovery of EC 22081— 1916 as a fast rotating, single sdB star of 
low gravity. Its atmospheric parameters indicate, that the hydrogen envelope 
must be unusually thick, which is at variance with the He- WD merger scenario, 
but consistent with a common envelope merger of a low-mass, possibly substellar 
object with a red-giant core. 

Subject headings: stars: individual(EC 22081 — 1916) — stars: horizontal-branch 



1. Introduction 

Hot subdwarf B stars (sdBs) are evolved core helium-burning stars wi th very thin hydro- 



gen envelopes residing at the extreme end of the blue horizontal branch (iHeberl |2009| ). The 
formation of these stars remains unclear. Normal horizontal branch stars are formed after 
core helium-burning is ignited in the red-giant phase. Since the hydrogen envelopes of sdBs 
are extraordinarily thin, large mass-loss is necessary right at the tip of the red-giant branch. 
In the case of the close binaries among the hot subdwa rfs - about half of the known sdB stars 



are members of short-period (P < 10 days) systems (IMaxted et al.ll200ll : iNapiwotzki et al. 
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2004al ) - the required mass loss is triggered by the formation of a common envelope, which 
is finally ejected. The formation of sdBs with main sequence co mpanions in wide orbi ts on 
the other hand can be explained by stable Roche lobe overflow (IHan et al.ll2002l . |2003[ ). 



The formation of single sd Bs is less wel l understood. Several single star scenarios are 
currently under discussion (see iHeben |2009| . for a review) , but all these scenarios require 
either a fine-tuning of parameters or extreme environmental conditions which are unlikely to 
be met for the bulk of the observed subd warfs in the field. A partic ularly interesting scenario 
has been suggested by ISokerl (Il998l ) and iNelemans fc TaurisI (119981 ) . Accordingly, single sdB 
stars could be formed if a substellar companion in close orbit was engulfed by the red-giant 
progenitor and provided sufficient angular momentum to the common envelope before it was 
destroyed. 

Alternative scenarios invoke stellar mergers to form single sd B stars. The me r ger o f 
binary white dwarfs was investigated by IWebbinkI (119841 ) as well as llben fc Tutukovl (119841 ) 
who showed that an EHB star can form when two helium core whit e dwarfs merge and the 
product is sufficiently massive to ignite helium. jPolitano et al.l ( l2008l ) proposed the merger of 
a red giant and a low mass main-sequence star during the common envelope phase. However, 
the merger chann el was under debat e, because all single sdBs analysed so far turned out to 
be slow rotators (iGeier et al.l l2009a] ) in contrast to the expectations. Here we report the 
discovery of the fast rotating, single sdB EC 22018— 1916. 



Observations 



EC 22018-1916 {V = 12.9 mag, aaooo = 22'^10'"52?9, 



2DQ0_ 



-19°01'50") was discovered 



in the course of t he Edinburgh-Cape blue o bject survey (IStobie et al.lll997( ) and classified 
as an sdB star by ICopperwheat et al.l (l201ll ). Five high resolution spectra were taken with 
the FEROS spectrograph [R = 48 000, A = 3750 -9200 A) mounted at the ESO/MPG2.2m 
telescope at La Silla. The first spectrum was taken on June 14, 2006 followed by two 
exposures on August 9 and 11 in the same year. The last two spectra were taken consecutively 
on October 30, 2010. In total, the data points cover a timespan of 4.5 years. The spectra 
have been reduced with the FEROS-DRS pipeline in the context of the MIDAS package. A 
median filter was applied to correct for cosmics. 

EC 22018— 1916 has been monitored by planetary transit surveys. Du e to its favourable 
declination it was observed by both the All Sky Au tomated Survey (AS A S, |Pojmanskilll997l ) 
and the Northern Sky Variability Survey (NSVS, IWozniak et al.ll2004l ). Both light curves 
have been downloaded from the data archives. The ASAS dataset contains more than 450 
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data points of Johnson-V photometry. We selected data extracted with smallest aperture 
and included only measurements of sufficient quality (flagged from A to C). The light curve 
was folded to a period of 0.15 d (see Sect. H]) and binned. The white light curve from NSVS 
contains 93 data points and was phased in the same way. No variations exceeding 2% (NSVS) 
and 1% (ASAS) were detected. 

The Vizier database contains several consistent proper motion measureme nts of this ob 



iect. Among them the PPMXL (/x„ cos 6 = 17.0 mas yr ^, fis = —13.5 mas yr ^.iRoeser et al. 
2010) and UCAC3 {fXaCosS = 17.7 mas yr~^, fxs = —12.1 mas yr~^, Zacharias et al. 201ol ) 



values are independently measured and perfectly consistent within the error bars. Taking 
the average values we obtain ^ = 21.6 mas yr~^. 



3. Atmospheric parameters and rotational broadening 



Atmospheric parameters and projected rotational velocity (see Tabled]) have been deter- 
mined simultaneously by fltting a grid of synthetic spectra, calculated from line-blanketed, 
solar-metalicity LTE model atmospheres (iHeber et al.ll2000l ). to the hydrogen Balmer lines 
(H;3-Hio) and helium lines (He 1 4 026,4472,4922,5876 A and He ll 4686 A) using the SPAS rou- 



tine developed by H. Hirsch (e.g. iGeier et al.ll2011bl ). The single spectra have been corrected 



for their orbital motion and coadded. Statistical errors are determined with a bootstrapping 
algorithm. 

As can be seen in Figures [1] and [21 the Balmer line cores and the helium lines are 
signiflcantly broadened. In order to flt these lines, a very high rotational broadening of 
f rot sin i = 163 ± 3kms^^ is necessary. The fact that no metal lines have been found is 
consistent with this high broadening since weak features melt into the continuum in this 
case. The resulting effective temperature Tes = 31100 ± 1000 K, and helium abundance 
log y = — 1.97±0.02 are typical for sdB stars, whereas the surface gravity of log g = 4.77±0.10 
is unusually low for the effective temperature in question (see Figure E]). 

The radial velocities of the flve single spectra were measured by fltting model spectra 
with flxed parameter s derived from the spectr al analysis to the Balmer and helium lines using 
the FITSB2 routine (INapiwotzki et al.ll2004bl . see Tabled]). No signiflcant RV variations were 
measured. The radial velocity of the star is constant at —13.1 ± 3.6kms~^. 

By adopting the canonical mass of sdB stars (0.47 Mq) we can derive the distan c e from 
the atmospheric parameters and the apparent magnitude following iRamspeck et al.l ( 120011 ). 
The transversal velocity Vt of the star in kms~^ can then be calculated using the simple 
formula Vt = A.YAdfi, where the distance d is given in kpc and the proper motion fi in mas. 
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The distance to the star is 1.5 kpc and the transversal velocity ^ 150 km s ^ per fectly 



consistent with an evolved star in the thick disk or in the halo (e.g. iTillich et al.ll201ll ). 



4. Constraining the nature of EC 22018-1916 

EC 22018— 1916 has the highest f rot sin i ever m easured for an sdB star. All other single 
sdB stars analysed so far have frot sini < lOkms"^ (iGeier et al.ll2009al ). In the following we 
discuss and exclude several possible explanations for this finding. 



Main sequence star? Rotational velocities exceeding 100 km s~^ are quite common among 
main sequence A and B stars. Since the surface gravity logg = 4.77 is at the lower end of the 
hot subdwarf parameter range, the star may be regarded as a missclassified massive main 
sequence star. This interpretation, however, can be ruled out because the surface gravity is 
too high (see Figure [3]) and the helium abundance (1/10— solar) far to low. A double-lined 
binary consisting of two hot main sequence stars is another option. An unresolved, double- 
lined binary may explain the high measured surface gravity, which could be overestimated in 
this case. However, we cannot imagine a combination of main sequence stars, which would 
produce such an unusual spectrum. 



Hot subdwarf with unresolved pulsations? Two kinds of sdB pulsators are known. 
The slow pulsations of the V 1093 Her stars (sdBVg) are not expected to influence the line 
broadening significantly. In the case of the short-period pulsators (V361Hya type, sdBVr), 
unresolved puls ations can severely affect the broadening of the lines and therefore mimic 
higher v^ot sin i. iTelting et al.l ( l2008l ) showed that this happens in the case of the hybrid pul- 
sator Balloon 090100001. Unresolved pulsations are also most likely responsible for the high 
39kms~^ measured for the high-amplitude pulsator PC 1605+072 ( iHeber et al. 



sm I 



19991). 



The five spectra of EC 22018— 1916 have been taken with exposure times ranging from 
900 s to 1500 s. The effective temperature of the sdB is consistent with the ones of short- 
period pulsating sdBs. The typical pulsation periods of sdBVr stars are of the order of a few 
minutes and therefore shorter than the exposure times. However, the measured f rot sin i = 
163kms~^ is so high that very large photometric variations at periods of a few minutes 
would be inevitable. Therefore the broadening cannot be caused by unresolved pulsations. 
The prominent mode of the strongest known sdB pulsator PC 1605+072 has a photometric 
amplitude of ~ 13% (IKoen et al.lll998[ ) and an RV amplitude of ~ 15kms~^ (lO 'Toole et al. 



- 5 - 



2005al ). In order to cause the line broadening necessary to fit EC 22018— 1916, both values 
would have to be much higher. Due to the fact that we neither detect RV variations nor 
any features in the light curves (see Figure S]) we conclude that EC 22018— 1916 is not a 
high-amplitude pulsator. 



Hot subdwarf with high magnetic fields? lO 'Toole et al.l (j2005bl ) discovered magnetic 
fields up to ~ 1.5 kC in a small sample of sdB stars. From the analysis of magnetic white 
dwarfs with field strengths in the MG range it is known that sm all Zeeman splitting can 
mimic a broadening of the spectral lines (see e.g.. lKiilebi et al.ll2009l ). Could EC 22018— 1916 
be the prototype of a new sdB class with very strong magnetic fields? Looking at Figs. [Hand 
[2] this explanation can be ruled out as well. Zeeman splitting affects every single spectral line 
in a different way. In contrast to that, the broadening of the lines is uniform. We therefore 
conclude that the line broadening of EC 22018— 1916 is caused by rotation. 



Close binary with large RV amplitude and orbital smearing? The hot subdwarfs 
with the highest measured projected rotational velocities (f rot sin i > lOOkms"^) all reside 
in very close binary systems with orbital periods of ~ 0.1 d. These sdBs were spun up by 
the tidal influen ce of their close compan ions and their rotation became synchronised to their 
orbital motion JCeier et al.lboOTi boioh . 



A close companion would therefore be the most natural explan ation for the high f rot s in i 
of EC 22018-1916. The colours of this s tar (J-Ks ~ 0.0, 2M ASS: Iskrutskie et "al1l2006[ ) do 
not show any signs of a cool companion (jStark &: Wadell2003[ ). A possible unseen companion 
must therefore be either a low mass main sequence star, a compact object like a white dwarf 
or a substellar object. 

White dwarf and main sequence companions can be immediately ruled out because no 
significant RV variations are detected on timescales of years, days and half an hour. This is 
also a strong argument against the hypothesis, that the strong line broadening may be at 
least partly caused by orbital smearing. Since the exposure times of the FEROS spectra are 
rather long (900 — 1500 s), the RV shift during t he exposure would be large in a close binary 



with high RV amplitude (e.g. iGeier et al.l 120071 ). However, no RV variations are measured 
(see Figure IH upper panel). Furthermore, a quantitative spectral analysis as outlined in 
Sect. |3]has been performed for the single FEROS spectra and no significant variations in the 
atmospheric parameters or the f rot sin i were detected (see Tabled]). Orbital smearing can 
therefore be ruled out. 
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Close binary with small RV amplitude? The remaining option would be a sub- 
stellar companion in a very close orbit similar to the sdB+BD binary SDSS J08205+0008 
( IGeier et al.ll2011al ). Assuming that the rotation of EC 22018— 1916 is synchronised an up- 
per limit for the orbital period can be calculated. Adopting MgdB = 0.47 Mq and using the 
measured \ogg the radius of the star is -RsdB = ^/ M^dsG/g — 0.47 i?©. Taking the incli- 
nation into account (frot > "^rotsini) we can calculate an upper limit for the orbital period 
P < 2TTRsAB/vrot sini — 0.145 d. Another strict constraint is set by the lack of significant RV 
variations. Taking the standard deviation of the RV measurements and multiplying it by 
three we end up with a conservative upper limit for the RV semiamplitude K < 12kms^^. 

Adopting the upper limits for P and the companion would have to be a brown dwarf 
with ^ 20 Mj and a radius of ^ 0.1 i?©. For inclinations lower than 90° the orbital period 
of the putative binary must be shorter, because the absolute rotational velocity of the sdB 
has to be higher to keep frot sini fixed at the observed value. 

However, other important constraints have to be met as well. Neither the sdB nor its 
putative companion are allowed to fill their Roche lobes, because in this case the system 
would exchange mass. Since no indicative features for ongoing mass transfer (e.g. emission 
lines) are present in the spectra, t he system must be detached. Calculating the Roche radii 
of both components as outlined in lEggletonl (Il983[ ) we derive a minimum orbital period for 
the system of ~ 0.11 d and a minimum inclination of 47°. For shorter periods and hence 
lower inclinations the sdB would fill its Roche lobe. A similar limit (~ 0.1 d) is derived, if 
we allow K to be smaller than 12kms~^. For orbital periods shorter than that, a putative 
brown dwarf companion would fill its Roche lobe. 

These simple calculations show that the possible parameter space of a close and synchro- 
nised binary would be extremely narrow (P ~ 0.1 — 0.15 d, i^' ~ 4 — 12 kms~^). Furthermore, 
all possible configurations would lead to photometric variabilities easily visible in the light 
curve. Close sdB+dM or BD systems are not only often eclipsing, but also show sinusoidal 



varia t ions due to light from the irradiate d surface of the cool companion (e.g. l0stensen et al. 



2010; 


For et al. 


2010; 


Geier et al. 


2011a) 



ity EC 22018— 1916 has a very high luminosity compared to other sdBs, which should lead 
to a very strong reflection effect at inclinations of 50° or higher. Since no variations were 
found in the ASAS and NSVS light curves (Figure H]), a nearby low- mass companion can be 
excluded as well. 
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Conclusion 



After excluding all possible alternative scenarios we conclude that EC 22018— 1916 is 
the first single sdB star which is rapidly rotating. Furthermo re, the logq of EC 2201 8— 1916 
is the lowest one ever measured for an sdB (see Figure [3]) . I0stensen et al.l (120 111 ) argued 
that the low gravity of the pulsating sdB J20163+0928 (log q = 5 .15) may be due to a rather 



thick layer of hydrogen. In the model of iHan et al.l (120021 . |2003| ) even the merger remnants 
with the highest masses would need a hydrogen layer of ~ 0.01 Mq to reach at such low 
surface gravities. 

The formation of such an object through single star evolution is very hard to explain. 
EC 22018—1916 thus might have been formed by a mer ger event. Three merg er scenarios have 
been p roposed to explain the origin of hot subdwarfs. IWebbinkI (Il984j ) and llben &: Tutukov 
(119841 ) proposed th e merger of two He-WDs as possib l e formation channel, which has been 
further explored by ISaio fc Jefferyl (120021 ) . iHan et al.l (120021 . |2003| ) included this channel in 
their bin ary evolution cal culations and were able to model both the UV-excess in elliptical 
galaxies (IHan et al.l 120071 ) and t he differen t close binary fractions of sdBs in populations of 
different age in a consistent way (lHanll2008l ). He- WD mergers are believed to have very small 
envelope masses and are expected to be situated at the very blue end of the EHB. Both is 
at variance with the p osition of EC 22018— 1916 in the (Tgfj — log (?)- diagram (see Figure [3]). 
Justham et al.l (120101 ) proposed that the merger of a close binary system consisting of an 
sdB and a He- WD may form a single helium enriched sdO. EC 22018— 1916, however, is 
helium-deficient . 

EC 22018— 1916 most likely belongs to an old stellar population, either thick disk or 
halo. Its position in the (Teg-log 5')-diagram (see Figure [3]) ma y indicate a mass highe r than 
canonical, which would be consistent with the predictions by iHan et al.l (120021 . 120031 ). If it 
should be the remnant of a He- WD merger, this would imply important constraints on the 
merger process itself. Since the helium abundance of EC 22018— 1916 is ten times below the 
solar value, enough hydrogen must have survived the merger and must have been enriched 
in the atmosphere by diffusion processes. 



T he thi rd channel was suggeste d by lSokerl (119981 ) and further explored by lSoker fc Harpaz 



( I2OOOI . I2OO7I ). iPolitano et al.l (l2008h followed this idea and focused on the formation of hot 
subdwarfs. The merger of a red-giant core and a low-mass, main-sequence star or substellar 
object during a common envelope phase may lead to the formation of a rapidly rotating hot 
subdwarf star. This scenario fits particularly well with observations for several reasons. 



First, the helium core of a red giant merges with an unevolved low- mass star or a brown 
dwarf. Both have hydrogen-rich envelopes. This provides a natural explanation for the low 
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He abundance and surface gravity of the remnant. The hydrogen is provided by the merged 
companion. Furthermore, this companion also provides the energy required to eject the 
envelope and form the sdB. Several sdBs with low mass stellar and substellar companions 
have been f ound most recently and the true number may be much higher due to selection 
effects (e.g. IPor et al.l[201o[ lOstensen etHlboid ICeier et al.lboOQbl . I2OI1/ 



A very important predicition made by lPolitano et al.l (120081 ) is that sdBs formed via the 
CE-merger channel should be rare. EC 22018 — 1916 is unique amon g ^ 100 slowly rotating 
sdB stars analysed so far (iGeier et al.ll2009al ). In contrast to that, iHan et al.l ( 12002| . 120031 ) 
predict a large fraction if not all of the single sdBs to be formed by WD-mergers. Unless there 
is a mechanism to get r id off all the angular momentum involved in a merger as suggested 
by lSaio fc Jefferyl (|2002[ ). this observation is hard to explain. 



Furthermore, iPolitano et al.l ( 120081 ) predict that a large fraction of the sdBs formed 
after CE-merger should rotate with a critical velocity fcrit, which is defined as the rotational 
velocity at which mass loss induced b y centrifugal forces pr events the red-giant core to accrete 
more material from the secondary. iPolitano et al.l (120081 ) estimate this critical velocity to 
be about one third of the breakup velocity fbr = (GM/i?)^/^. Using the parameters derived 
for EC 22018— 1916, we calculate fcrit — 145kms~^ perfectly consistent with the projected 
rotational velocity measured from the spectrum. 



In conclusion, the scenario proposed by ISokerl ( 119981 ) and IPolitano et al.l (120081 ) fits best 
with the observational data obtained so far, although the He-WD-|-He-WD or sdB+He-WD 
merger scenarios cannot be ruled out. EC 22018— 1916 is the first candidate for a merger 
remnant among the hot subdwarf stars. Similar objects are expected to be found in large 
spectroscopic databases like SDSS. Due to the high rotational broadening the quality of these 
data should be sufficient to find them. 



Based on observations at the La Silla Observatory of the European Southern Observa- 
tory for programmes number 082.D-0649 and 084.D-0348. S. G. is supported by the Deutsche 
Forschungsgemeinschaft (DFG) through grant HE1356/49-1. We thank L. Morales- Rueda 
for sharing her data with us. Furthermore, S. G. wants to thank Ph. Podsiadlowski, C. S. 
Jeffery, R. H. Ostensen and S. J. O'Toole for defending the merger channel as possible for- 
mation scenario for hot subdwarfs. Special thanks go to the organizers of the 4th sdOB 
meeting in Shanghai where these and other problems were discussed. 
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Table 1: Parameters of EC 22081-1916 
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Spectrum with low S/N. 
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Fig. 1. — Fit of synthetic LTE models to some hydrogen Balmer hnes. The thin sohd 
hne marks models without rotational broadening, the thick solid line the best fitting model 
spectrum with Vrotsini — 163kms~^. 
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Fig. 2. — Fit of synthetic LTE models to helium lines (see Figured]). The extreme rotational 
broadening of the lines is obvious. 
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Fig. 3. — Teff — log g diagram. The grey circles mark sdBs from the SPY proj ect (ILisker et al. 



2005h. The low gravity sdBs HS 1736+8001 (lEdelmann et all [20031), PG 1605+072 
jHeber et al.lllQQflh and J20136+0928 jOstensen et al.ll201ll ) are plotted as triangles. The 
helium main sequence (HeMS) and the EHB band (limited by the zero-age EHB, ZAEHB, 
and the terminal-age EHB, TAEHB] are superimposed with EHB evolutionary tracks for 
solar metallicity from lDorman et al.l (11993 ). The location of the main sequence (MS, limited 
by the zero-age MS, ZAMS, an the terminal-age MS, TAMS) is indicated by the long-dashed 
horizontal lines. 
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Fig. 4. — Time resolved photometry (AS AS V-band, NSVS white hght, lower panels), ra- 
dial and projected rotational velocities (with formal la fitting errors, upper panels) plotted 
against phase. The data has been folded to the upper limit for the orbital period of a puta- 
tive close binary with substellar companion (P ~ 0.15 d, see Sect. 4). Two complete phases 
are plotted for better visualisation. Since no significant variations can be seen in the data, 
EC 22018— 1916 is neither a strong pulsator nor in a close binary system. 



